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ABSTRACT 

We investigate rapid mass supply process by nested bars in the Galaxy by 
numerical simulation. We simulate gas flow in the whole galaxy disk with nested 
bars, which are the outer bar and the inner bar, especially with highly spatial 
resolution in the galactic central region. We assume two cases of inner bar size 
which are a smaller one and a larger one than the radius of the 200 pc gas ring 
which is corresponds to the Central Molecular Zone. From our numerical results, 
in the large size bar cases, the inner bars with large elongation induce sufficient 
mass inflow and destroy the 200 pc gas ring. On the other hand, in the small 
size bar cases, the inner bars with large elongation induce large mass inflow and 
do not destroy the 200 pc gas ring. This mass inflow is caused by straight shocks 
excited by the inner bar. In this case, nuclear gas disks of ~ 15 pc radius are 
formed. The nuclear gas disks are self-gravitationally unstable and we expect 
formation of compact star clusters under strong tidal force in the nuclear gas 
disks. We discuss evolution of the nuclear gas disk. 

Subject headings: Galaxy: center — Galaxy: kinematics and dynamics — meth- 
ods: numerical 



INTRODUCTION 



Gas fueling to a galactic center is very important for activity of active galactic nuclei, 
growth of supermassive black holes (SMBHs), nuclear starbursts, formation of super star 
clusters in a galactic central region, and other interesting phenomena. Our galaxy is very 
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interesting for this study because of the following reasons. First, the center of our galax y 
is the most closest galactic center. Its distance is about 8.0 kpc (jEisenhauer et al.ll2003l ). 
Therefore, there are many observational data with high resolution over wide wavelength. It 
is easy to compare those with numerical simulations. Secondly, there is evidence of recent 
mass supply in our galactic center. There are young massive compact star clusters (the 
Arches, Quintuplet, and Central clusters) in the Galactic center. The se clusters a re located 



within 30 pc from the Galactic center and have a number of OB stars (lFigerll2002l ). Because 



of the short age of the young stars, such massive st ar formation occurr ed within last several 



million years in the central region of the Galaxy (IMezger et al.lll996l ). Formation of these 
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(10 5 cm -3 ), clumpy, and turbulent with large li ne widths (> ^kms" 1 ) 
has radius of a few pc and a mass of « 10 6 M Q fjChristopher et all 120051 ) 
found a gas stream from the giant molecular cloud (the 20kms _1 cloud) near the Galactic 
center to the CND. This may be gas inflow to the CND. 

It is expected that a vast of gas is supplied from the central molecular zone (CMZ), which 
is a ring-like gas distributio n and extends over the range of Galactic longitud e —1.5° < / < 



19961). The size of 



2°, to the Galactic center (jSerabyn fc Morris! 119951 ; iMorris fc Serabyn 
the CMZ is ~ 200 pc. It is certainly form ed by the large-scale bar (jBinney et al.l 11991 
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1996 



Sawada et al.l 120041 ) and has a large amount of molecular gas 5- 



10 x 10 7 M g ( jSerabyn fc Morris 19951 ). However, it is unclear how the gas is transported 



further in . Secondary effect s like dissipation can then drive the gas fur ther in but at slowe r 
rate (e.g. Heller et al. 200lh. or likewise gravitat ional instabilities (e.g. Fukuda et al. 2000 ). 
magnetic viscosity (e.g. IMorris fc Serabynlll996l ). 



For the gas feeding, many authors show an important role of a bar (e.g. lAthanassoula 



19921 ). Nested bars, which consists of a outer bar and inner bars, may play an important role 
in the gas feeding to galactic centers. This idea was firstly proposed bvlShlosman et al.l (Il989l ) 
as a mechanism for fueling AGNs. I nspired by the idea of IShlosman et al. I dl989h. many nu- 
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Shen fc Dibattistall2007l ). iFriedli fc Martinetl (119931 ) performed three dimensional sim- 
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ulations of gas and stars and showed that an inner bar can drive gas infall to a galactic 
center. 



Nested bars are observed in nearby barred galaxies in a large fraction (~ 30%) (jWozniak et al. 



19951 : IFriedli et al.lll996l : lJungwiert et all 199 71 : lErwin fc Sparkdl2002l : lErwinll2004f ) . The large 
fraction indicates that nested bars are dynamically stable or recurrent structures. Nested 
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bars are expected to be dynam ically decoupling, since the orientations of both bars are 
random (IButa fc Crocker! 11993). Dynamical decoupling of these was also reported in many 
numerical studies (Friedli fc Martinet 19931; Macieiewski & , Sparke 2000 ; Rautiainen et al. 



20021 ; lEnglmaier fc Shlosman 



2004 



Debattista fc Shenl 120071 ) 



An increasing number of observational stud ies show effect of nested bars in gas flows 
in central regions of galaxies. iFathi et al.l (120061 ) observed the central region of the double- 
barred galaxy, NGC 1097, with high resolution, using GMOS-IFU and HST-ACS. They show 
clear evidence of radial streaming motion down t o abo ut 10 pc from the nucleus by mapping 
the gas velocity fields. ISchinnerer et al.l (120061 . 120071 ) observed molecular emissions in the 
central region of the nearby double-barred spiral galaxy NGC 6946 with very high spacial 
resolution (< 1") with the IRAM Plateau de Bure interferometer. They showed that there are 
nuclear massive gas clumps and straight dust lanes inside the inner bar. They co ncluded that 



the in ner bar is closely related with the pile-up of molecular gas to the nucleus. iMeier et al. 



(120081 ) observed the central region of barred galaxy, Maffei 2, with BIMA and OVRO. They 
show a nuclear ring, whose radius is ~ 80 pc and mass is 6.9 x 10 6 M Q , well inside the bar 
and that overall morphology of gas, including the nuclear ring, can be explained by a nuclear 
bar by comparing the position-velocity diagram of molecular gas with orbits of molecular 
clouds in their nuclear bar model. These studies support the important role of inner bars in 
transporting gas to a galactic center. 

Recently, observational studies show evidence of an inne r bar in our galaxy, which is 
much smaller than the outer bar of semi-major axis 3.5 kpc. lAlardl (120011 ) studied surface 
density of old stellar pop ulation in the inner bulge by using the 2MASS data and show 
evidence of the inner bar. iNishiyama et al.l (120051 . 120061 ) investigated the shift of the peak 
position of red clump stars distribution over \l\ < 10.5° using the IRAF 1.4 m telescope with 
the near-infrared camera SIRIUS and showed that the gradient of this shift clearly changes 
in \l\ < 4°. They interpreted that this structure may be due to the inner bar. 

We study the possibility that the inner bar play an important role in the mass supply 
from the CMZ to the Galactic center. Previous theoretical studies have no t reported the case 



of large contribution of inner bars in gas supply to a galactic center (e.g. iMaciejewski et al. 
20021 ; iRautiainen et al.l 120021 ). They studied limited cases. We investigate various inner 
bar models in this paper. We perform two dimensional hydrodynamical simulations in a 
gravitational potential model of our galaxy, assuming several inner bar models. In the 
simulations, we systematically change the mass and the axial ratio of inner bar models, since 
parameters of the inner bar is not clear from observations. 

In §2, we give our gravitational models and numerical method. In §3, we show the 
results of our simulations. In §4, we discuss gravitational instability and evolution of nuclear 
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gas disks, which are obtained in our numerical results. In §5, we summarize our study. 



MODEL 



2.1. Gravitational Potential of the Galaxy 



As t he gravitationa l pote ntial of our galaxy except for an inner bar, we assume the 
model of iBissantz et al.l (120031 1 for the Galactic b ulge, the stellar d i sk, th e outer bar, the 
spiral arms, and the dark halo (R > 500 pc) and lLaunhardt et al.l (120021 ) for the nuclear 
bulge (R < 500 pc) and the SMBH. lBissantz et al.l (120031 ) simulated gas motion in our galaxy 
potential model, which consists of the Galactic bulge, the stellar disk, the outer bar, the spiral 
arms and the dark halo. They gave pattern speed of the outer bar and the spiral arms (floB ~ 
60 Gyr^ 1 and f2sp ~ 20Gyr -1 , respective ly) to reproduce observational gas kinematics of 



molecular clouds. lLaunhardt et al.l (120021 ) analyzed IRAS and COBE DIRBE data of the 
central 500 pc of our galaxy. They gave mass distribution of the nuclear bulge, which is 
distinguished from the Galactic bulge by its flat disk-like feature, assuming a constant mass- 
to-light ratio. They estimated that the nuclear bulge has a mass of 1.4 ± 0.6 x 1O 9 M . We 
assume the rotation curve obtained from the mass distribution of the nuclear bulge and the 
SMBH in R < 500 pc for the rotation curve of the total stellar mass. We connect smoothly 
the rotation curves obtaine d from the nuc l ear bu lge and the SMBH in R < 500 pc and from 
the stellar component of of IBissantz et al.l (120031 ) in R > 500 pc. Details on th e gravitational 



poten tial of the outer bar, the spiral arms, and the dark halo are described in IBissantz et al. 



(120031 ). Fig. [U shows the rotation curve of one of our models, the model S33. In this figure, 
we use axially averaged mass distribution of the inner and outer bars. In Fig. [2l we show 
the angular velocity curve of the model S33. In this figure, there is the local maximum of 
Q — k/2 at 150 pc, where Q is the angular velocity and n is the epicyclic frequency. We 
point out that the curve of f2 — k/2 in R < 500 pc is rather uncertain, since it is difficult to 
measure accurately the mass profile in this scale. 



2.2. Inner bar potential 

We assume Ferrers bar models for inner bars, since a density profile of the inner bar is 
not observationally confirmed. The Ferrers bar model has a density distribution as 

/ ^2 

p(x,y,z) = Po (1-^--^-^ 
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where po is density at the origin (jFerrerslll877l ). We assume n — 1 and b = c. po is related 
with mass of the inner bar Mrs through p = g 5J ^f for n — 1. Parameters we choose are 
given in Table{T]and[2J 



We assume two cases of the length of the semi - major axis of t he inn er bar models, 
aiB = 200 pc and 600 pc from the following studies. IWozniak et al.l (119951 ) performed the 
BVRI survey of 36 disk galaxies selected as candidates for having an inner bar or a triaxial 
bulge within the outer bar. They showed that outer to i nner bar axis r atios, oob/ciib, are 
in the range of 3.7 to 18.0 with the mean value of 7.2. iFriedli et al.l (119961 ) observed 13 
disk galaxies, which had been classified into galaxie s likely having an inner bar or a triaxial 
bulge within the outer bar in IWozniak et al.l (119951 ). with JHK band. They show a similar 
result, 4.0 < aoB/ a iB < 13.4 with the mean value of 7.2. In our galaxy, the ranges above 
corresponds to oib = 200-875 pc for the semi-major axis of the outer bar, 3.5 kpc. If our 
galaxy is a normal nested barred galaxy, our assumed values of are is in this range. 

In our assumption on sizes of t he inner bar models, we also consider the fact that inner 
bars of ten coexist with n u clear rings (IButa &: Crockerlll993uShaw et al.lll993uErwin &: Sparke 
20021 ). lErwin fc Sparkei (120021 ) found that 60% of their sample galaxies with nuclear rings 
have inner bars. In such galaxies, inner bars are often surrounded by nuclear rings and the 
size of the inner bars is comparable with that of the nuclear rings. In our galaxy, if the CMZ 
corresponds to such a nuclear ring, the size of the inner bar may be comparable with the 
size of the C MZ (R « 200 pc). This is consistent with the projected size of the inner bar of 
Alardl (bpOlh that is ~ 1.5° — 2° « 200 — 300 pc. However, the size of the inner bar proposed 
by lNishiyama et al.l (120051 . 120061 ) is ~ 520 pc, and it is much larger than the size of the CMZ. 



Our ass umption on sizes o f the i nner bar models is con s istent with recent numerical 
simulations. iDebattista fc Shenl (120071 ) and lShen &: Dibattistal (120071 ) investigated formation 
of long-lived inner bar from a psuedobulge by performing N-body simulations. They showed 
that inner bar e nds are much smal l er tha n their corotation radius -Rcr- Similar result is 
also obtained in IFriedli fc Martinetl (ll993l ). The Rcr of the inner bar is as large as 600 pc 
in our models, if the pattern speed of the inner bar is near the local maximum of Q — k/2 
(our choice is intended to be consistent with the N-body simulations; see below). <zib may 
be less than 600 pc. Since the curve Q — k/2 in R < 500 pc is rather uncertain, we assume 
two cases of oib- Hereafter, we call the inner bar models with am = 200 pc small inner bars 
and with ctrB = 600 pc large inner bars. 

We assume that th e mass of the inne r bar ( see Sec. 12.31) is a part of the mass distribution 
of the nuclear bulge of lLaunhardt et al.l (120021 ). We give the mass of the inner bar models 
in Table. [1] and [2j As shown in Fig. [Tj, the mass of the inner bar models are quite smaller 
than the total mass within the radius of semi-major axis of the inner bars. 
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We assume that pattern speeds of the inner bar models are near the local maximum 
of Q — k/2, which is located at about 150 pc (se e Fig. [2]). This is consistent with N-body 



simul ations of formation of nested barred galaxies fjFriedli fc Martinetlll993l ; lRautiainen et al. 



20021 ). We also assume that the inner bars are prograde. In some small inner bar models, we 
change the pattern speed around the local maximum of Q — k/2 to investigate the effect of 
the pattern speed on mass inflow rate to the galactic center. The range of the pattern speed 
in each model is 175-375 km s _1 kpc -1 . We give the pattern speeds in Table [3j 

We use Qt = {F(t>j -frjmax as a measure of the strength of the inner bar, where {F^/ F r ) rDa _ x 
is a maximum value of a ratio of the azimuthal component of the gravitational force of the 
inner bar to the radial component of the total gravitational force within 500 pc. Qt of each 
model are given in Table [TJ and [2J 

Hereafter, we specify the inner bar models by the name in the Table [TJ and [2] and by the 
value of the pattern speed as S33 (Qib = 300 km s" 1 kpc -1 ). In the case of the large inner 
bar models, we omit the value of the pattern speed, since we assume the same pattern speed 
for them. 



2.3. Numerical method 



We use the advec tion upstream splitting method (AUSM) for numerical hydrodynamics 
(ILiou fc Steffenl Il993l ) . The AUSM is one of flux vector splitting schemes. In the AUSM, 
advection and propagation of acoustic wave are recognized as physically distinct processes. 
Therefore, the advective terms and the pressure terms in the flux vector are split separately. 
This makes the formula of the flux vector at the cell face very simple and leads to a reduction 
of numerical operations without loss of accuracy. The robustness and good performance of 
the AUSM in the app l ication to galactic gas s i mulations are we ll tested by many authors 
JColina fc Wadal I2OO0I ; IWada fc Normanl l200ll ; iMori et al.l I2OO2I 1 To obtain higher order 
spatial resolu tion, we use the second o rder MUSCL interpolation with the van Albada limiter 
function (e.g. iRadespiel fc Krolllll995l ). The AUSM with the MUSCL interpolation is easy to 



implement due to its simple form and well suitable to capture shock w aves in even rarefied 



medium. In our simulations , we do not use "gas recycling law" (e.g. lAthanassoulal 11992 



Englmaier fc Shlosmanl|2000l ). since we do not intend to seek a steady state of the flow and 
our simulation time is much shorter than a timescale of exhaust a large fraction of gas in 
the systems by star formation. 



In order to resolve gas motion in the galactic center region, we use two dimensional polar 
grids extending from 5 pc to 10 kpc in the Galactic radius. We divided radial grids into 370 
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logarithmically and azimuthal grids into 300 equally keeping the shape of each cell nearly 
square. The radial spacing AR of the grids decreases inwards. Very high spacial resolution 
is achieved in the central region, e.g. AR « 0.1 pc at R = 5 pc. 

We assume isothermal, non-self-gravitating, and non-viscous gas for simplicity. We 
do not consider a viscous term in the hydrodynamical equations. We use the equation of 
state of ideal gas with temperature of 10000 K, which corresponds to the sound speed of 
c s ~ lOkms" 1 and random motion of interstellar gas implicitly. We do not consider star 
formation and feedback process, such as supernovae and stellar mass loss, in this paper. 

We assume a rotationally supported gas disk for the initial state. This disk is flat and 
has infinitesimal thickness. Its outer radius and mass are lOkpc and 1O 1O M , respectively. 
The initial surface density of the disk is uniform in all models. 

The radial outer and inner boundary conditions are free and the azimuthal boundary 
condition is periodic. We record mass flux passing through the inner and the outer boundary 
for checking mass conservation. 

In order to avoid spurious phenomena, we introduce the non-axisymmetric components, 
such as the inner bar, the outer bar, and the spiral arms, of the gravitational potential slowly, 
compared to the rotational speed of each component. We gradually deform the gravitational 
potential of the inner bar from a spherical shape, 

(X 2 -|- y 2 -\- z 2 \ 
1 r§ / ' ^ 

where p' n = ^r^pr a nd r = to its assumed one from t = 100 Myr to 250 Myr as in 



Athanassoulal (119921 ). We also similarly introduce the Fourier component of t he gravitational 



potential of the outer bar and the spiral arms given by lBissantz et al.l (120031 ) from t = Myr 
to t = 100 Myr. 

We use the super computer SR11000/K1 of the Hokkaido university Information Initia- 
tive Center (IIC) for our simulations. 



3. NUMERICAL RESULTS 

We perform the hydrodynamical simulations for various masses and axial ratios of the 
inner bar systematically. In the small inner bar models, we also vary their pattern speed. 



We find a large amount of gas concentration to the galactic center in both sizes of the 
inner bar models. In small inner bar models, the inner bars induce gas inflow to the galactic 
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center for 0.05 < Qt < 0.3 without destroying the 200 pc gas ring, if Qjb ~ (f2 — K/2) max ~ 
SOOkms^kpc" 1 for 0.12 > Q T > 0.3 and if Q 1B w 225 km s" 1 kpc -1 for Q T > 0.05. On the 
other hand, in large inner bar models, gas concentration occurs if Qt > 0.1 and the inner 
bar destroys the 200 pc gas ring. In the following subsections, we describe the results in 
more detail. 



3.1. The no-inner bar case 



We perform hydrodynamical simulation in the Galaxy model without the inner bar to 
compare with the models with the inner bars. We show result of the no-inner bar model 
(model N) in Fig. [31 

Gas ridges are formed in the outer bar region by t = 100 Myr. Gas in the galactic disk 
flows into the central reg i on (R < 300 pc) along the gas ridges. This result is almost the same 



results of iBissantz et al.l (j2003l ). In our numerical results, a gas ring is formed at the radius 
of 150-250 pc. Its mass is almost constant at the value of ~ 3 x 10 8 M Q after t = 100 Myr. 
The mass and size of the ring correspond to the CMZ, of which extent is —1.5° < / < 2° and 
mass is 5-10 x 1O 7 M . We find similar gas rings in the models with inner bars. Hereafter 
we call these rings the 200 pc gas rings. 

The radius of the 200 pc gas ring is well inside the positi on of the ILR of the outer 
bar (~ 750 pc). This result agrees with iRegan fc Teubenl (120031 ). They showed that size of 
nuclear ring is related to population of X2 orbits, rather than the position of ILRs of an outer 
bar when gas motion is in the non-linear regime of hydrodynamics in the barred potential. 

Inside of the 200 pc gas ring, there are weak gas spirals. Their pattern speed agrees with 

the pa tt ern speed of the outer ba r . The se spirals are density waves found by lEnglmaier fc Shlosman 
( 120001 ). lEnglmaier fc Shlosman! (120001 ) show that gaseous spirals are formed inside the ILR 
of a bar in their numerical simulations of non-self-gravitating gaseous disks and that such 
gaseous spirals are supported by pressure force and stationary in the bar frame. The gaseous 
spirals in our simulation have similar property. Hereafter we call these spirals the nuclear 
spirals. The nuclear spirals become more tightly wound as approaching to the galactic cen- 
ter. Near 20 pc from the center, nuclear spirals are highly tight winding. An average mass 
inflow rate from 100 Myr to 500 Myr is very small, pa 3.6 x lO~ 4 M yr _1 . This radial gas 
inflow may be due to the nuclear spirals, since the total gravitational torque on the gas in 
the nuclear spirals region is consistent with the average mass inflow rate. In order to confirm 
this, we calculate the total gravitational torque on the gas inside R = 60 pc from the outer 
bar. The time averaged total gravitational torque within R = 60 pc between 100-500 Myr 
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is —7.6 x 10 49 g cm 2 s 2 . The mass inflow rate is as large as ~ 10 3 M yr 1 by this torque. 
This is consistent with the average mass inflow rate from 100 Myr to 500 Myr. 



Such nuclear spirals were not formed in the simulations of iBissantz et al.l (120031 ). This 
m ay be due to the lack o f the spacial resolution in the nu clear spirals region in the simulations 
of IBissantz et al.l (120031 ) . lEnglmaier fc Shlosmanl (120001 ) have shown that in simulations with 
insufficient spacial resolution to resolve nuclear spiral waves, they are quickly damped out 
due to numerical viscosity. 



3.2. The small inner bar models 

We find that a large amount of gas concentrates to the galactic center in the small inner 
bar models with Qt > 0.05 in some range of Vim- We divide our results into two cases, the 
high gas mass concentration case and the low gas mass concentration case. 

If Qt ^ 0.12 (S42, S43, S33, S34), high gas mass concentration to the galactic center 
occurs for both QrB ~ (fi — ft/2) max ~ 300 km s -1 kpc" 1 and fipB ~ 225 km s" 1 kpc -1 . If 
0.05 < Qt < 0.12 (S41, S32, and S23), high gas mass concentration to the galactic center 
occurs only for f2jB ~ 225 km s -1 kpc" 1 . One exception is the model S24, in which high gas 
mass concentration occur for both ~ (f2 — ft/2) max and f2iB ~ 225kms _1 kpc -1 in spite 
of 0.05 <Q T < 0.12. 



3.2.1. The high gas mass concentration cases 

Here, we describe the results of the model S33 (f2n3 = 300kms _1 kpc -1 ) in detail, since 
time evolution of gas distribution in the inner bar region are similar to the high gas mass 
concentration cases. 

We show the time evolution of the surface density of gas in the central 1 kpc square in 
the model S33 (f2n3 = 300 km s" 1 kpc -1 ) in Fig. HI One of characteristic gas distribution 
is straight shocks inside the inner bar (see Fig|5]). These shocks appear after the inner bar 
potential is introduced and become stronger with calculation time (see Fig. HH-f). These 
shocks extend from the galactic central disk to the inner edge of the 200 pc gas ring and are 
efficient to supply a large amount of gas to the galactic center. A massive nuclear gas disk 
forms in R < 15 pc. Its mass reaches as large as 10 7 M at t = 500 Myr. Hereafter we call 
this disk the nuclear gas disk. 

An elliptical gas ring is formed around the inner bar and is elongated along the inner 
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bar (see Fig. Hf). Similar elliptical gas ring is shown in iMaciejewski et al.l (120021 ) . Shape 
and surface density of this ring changes as the inner bar rotates. In Figj5l the ellipticity of 
the ring is larger at AO = 90° than at AO = 0°, while surface density of the ring is higher 
at AO = 0° than at A9 = 90°, where A9 is the angle between major axes of the inner bar 
and the outer bar. The velocity fields in the elliptical ring are smoothly connected to that 
of surrounding gas. 

In Fig. [7] we show the time evolution of the gas mass within 20 pc, M 2 o(t), in the 
model S33 (Qm = 300 km s -1 kpc -1 ). As deformation of the inner bar proceeds (t = 100- 
250 Myr), M2o(t) rapidly increases wit h time. Then, M^n i t) sat ur ates (t = 250 — 350 Myr) . 



Similar phenomenon was reported by IMaciejewski et al.l (120021 ). IMaciejewski et al.l (120021 ) 
showed that an inner bar keeps gas away from the galactic center and gas inflow due to the 
inner bar is negligible after it reaches its full strength. In the corresponding stage, in our 
simulations, velocity fields and gas distribution inside the 200 gas pc ring is perturbed by 
the inner bar. After velocity fields and gas distribution is quasi-steady, increase of gas inflow 
to the galactic center begins at t = 350 Myr. This second inflow continues to the end of the 
simulations. M2o(t) attains ~ 10 M Q at t = 500 Myr. An average second mass inflow rate 
is ~ 10 7 M n /100 Myr & 0.1 M yr _1 . We discuss the difference between our results and that 
of IMaciejewski et all (bo02h in Sect. 4. 



Occurrence of the second mass inflow depends on pattern speed of the inner bar. We 
show time evolution of M 2 o{t) for various pattern speeds of the inner bar in the model 
S33 in the lower panel of FigJTl Figure [7] shows that the second mass inflow occurs when 
the pattern speed is in 290-325 km s~* kpc" 1 and in 200-225 km s" 1 kpc" 1 in the model S33. 
We summarize M2o{t = 500 Myr) for small inner bar models in Table HI In this table, 
we denote models in which the second mass inflow occurs by bold letters. We note the 
models by daggers, in which M 20 (t = 500 Myr) exceeds the stellar mass within 20 pc M*(< 
20 pc) « 2 x 10 7 M (see Fig. 14 in lLaunhardt et al.l 120021 ) . In this case, we should solve 
self-consistently both the motion of gas and stars in the Galactic central region. Double 
daggers show the models in which the second mass inflows begin just before the end of the 
simulation. In these models, more mass will inflow into the galactic center, if we continue 
the simulations. 

Mass of the nuclear gas disk increases with time periodical l y by th e second mass inflow. 
Similar periodicity have been reported in IShlosman fc Hellerl (j2002al ). This case may be 
closely related with resonance phenomena between the outer bar and the inner bar. In our 
results, sufficient elongation for a small inner bar and a suitable Qib are needed for the 
second mass inflow. 
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3.2.2. The low gas mass concentration cases 

In small inner bar models with Qt < 0.05, gas mass concentration to the galactic center 
is small (S31, S21, S22, S13, and S14). In these models, loose two gas spirals appear in the 
inner bar region instead of straight shocks. These nuclear gas spirals become tightly wound 
near the center and a less massive gas disk appears in R < 20 pc from the center. 

The low gas mass concentration is due to the absence of the second mass inflow. We 
show the time evolution of M 2 o(t) in the model S21 = 300 km s" 1 kpc -1 ) by a dashed 
line in the upper left panel of Fig. [7J This figure shows that M 20 (t) saturates after the first 
mass increase and the second mass inflow does not occur till the end of the simulation. We 
test the time evolution of M2o(t) of the model S21 for various pattern speed, as shown in 
the upper right panel of Fig. [7J There is no second mass inflow in a range of Qib = 175- 
325kms _1 kpc -1 . Thus, we conclude that the second mass inflow needs Qt > 0.05. 

We address characteristic gas distribution seen in the low gas mass concentration models, 
since it is clear evidence of an inner weak bar. In Fig. [6l we show the snapshots of surface 
density of the model S21 for two pattern speeds, ilm = 200 and 300 km s" 1 kpc -1 . In both 
models, the loose two gas spirals are formed in the inner bar region and are surrounded 
by the gas rings. Similar structure is observed in the double-barred gal axy, NGC 1097. In 



this galaxy the l oose gas spiral s are observed with in the starburst ring (jPrieto et al.l 12005 



Fathi et al.ll2006l ). Contrary to iPrieto et all (120051 ) . it is possible that loose gas spirals are 



formed by an inner bar without any peculiar assumption. 



3.3. The large inner bar models 

3.3.1. The high gas mass concentration cases 

In large inner bar models with Qt > 0.11 (L42, L43, L33, L34, and L35), a large amount 
of gas concentrates to the galactic center for f2iB = 325 km s -1 kpc -1 . In Fig. [H] we show 
the time evolution of the surface density of gas of the model L42. In t = 100-250 Myr, 
elongation of the 200 pc gas ring increases. As can be seen in Fig. [8b-d, the 200 pc gas ring 
is highly elongated in t = 150-300 Myr. At t = 350 Myr, the 200 pc gas ring shrinks to less 
than R = 150 pc (see Fig. [Sf). Then, large part of the gas of the 200 pc gas ring rapidly 
concentrates into the galactic center and a very massive gas disk is formed at the center. 
The mass of the disk highly exceeds 10 8 M . The final value of M2o(t) is unreal because of 
the same reason described in Sect. 13.2.11 
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3.3.2. The low gas mass concentration cases 

In the case of Q T < 0.11 (L41, L31, L32, L22, L23, L24, L13, L14, and L15), a large 
amount of gas do not concentrate to the galactic center. The inner bar changes the shape 
of the 200 pc gas ring into more elliptical. The orientation of the deformed 200 pc gas ring 
is almost parallel to the inner bar. In these models, there is no enhancement of the mass 
inflow rate to the center. The average mass inflow rate over the simulation time is as small 
as the no-inner bar case. 



DISCUSSION 



4.1. Mass supply due to nested bars 

We have shown that mass supply process due to the nested bars is very efficient process 
by the numerical simu l ations . There are possible scenarios of the mass supply to the Galactic 
center. lAthanassoulal (119921 ) showed that gas ridges can reach a galactic center if a large- 
scale bar is yery strong. However, the ax ial ratio of the outer bar of our galaxy is ~ 3 



( jStanek et al.lll997l ; iRattenbury et al.l 120071 ). Hence, it is unlikely that ma ss supply to the 
Galactic center is due to the 'past' strong outer bar. iFukuda et al.l (120001 ) simulated self- 
gravitational instability of a nuclear gas ring and showed that a part of gas in the ring falls 
into a galactic center, since the gas transfers its angular momentum to a very massive clump, 
which is formed due to the fragmentation of the gas ring and subsequent mass accretion by 
surrounding gas. This process can explain mass supply to the galactic center if the CMZ 
corresponds to such a nuclear gas ring. In this simulation, as the result of the fragmentation, 
the nuclear ring is disrupted. This is not consistent with the CMZ in our galaxy. 

We have shown a large amount of gas concentration to the Galactic center, by performing 
two dimensional hydrodynamical simulations with various inner bar parameters (size of semi- 
major axis, mass, axial ratio, and pattern speed of the inner bar). We have performed 
simulations for inner bars with fi IB ~ (fl — K/2) mFUC , since th i s pattern speed is cons istent 



with the N-body simulation results (IFriedli fc Martinetlll993l ; iRautiainen et al.ll2002l ). We 
also have performed simulations changing the pattern speed of the inner bar for the small 
inner bar models to investigate effect of the pattern speed on mass inflow rate. We have 
assumed the two sizes of the semi-major axis of the inner bar, 200 pc and 600 pc. We have 
found the high gas mass concentration in both size of the inner bar. 

In the small inner bar models, The high gas mass concentration occurs for certain ranges 
of Qt and Qib- In the models with Qt > 0.12, the second mass inflow to the galactic center 
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occurs for Q\b ~ (f2 — ft/2) max . However, in models with Qt < 0.12, the second mass inflow 
does not occur for £7ib ~ (f2 — ft/2) max . For 0.05 < Qt < 0.12, the second mass inflow occurs 
for f2jB ~ 225 km s _1 kpc" 1 . Thus, the high gas mass concentration cases for the small inner 
bar models are divided into two cases: 

1. 0.05 < Q T < 0.12 and Q m ~ 225 km s" 1 kpc" 1 

2. Qt > 0.12, and fi ffi « (0 - re/2) max or fi IB ~ 225 km s" 1 kpc" 1 

One exception is the model S24, in which high gas mass concentration occurs for both 
n lB « (fi - K/2) max and fi ffi ~ 225 km s" 1 kpc" 1 in spite of 0.05 < Q T < 0.12. These 
results are summarized in Fig. [9j In this figure, the results of the model S25(f2iB = 
250 km s" 1 kpc" 1 ) and the model S41(f2iB = 250 km s" 1 kpc" 1 ) occupy the same point at 
(Qt, ^ib) = (0.115,250). High gas mass concentration occurs in the model S25(f2iB = 
250 km s" 1 kpc" 1 ), while it does not occur in the model S41(firB = 250 km s" 1 kpc" 1 ). 

The second mass inflow rates change periodically in the models which are denoted by 
asterisks in Table H] (see also the lower panel of Fig. [7]). These periodic changes imply that 
the second mass inflow is a resonance phenomenon between the outer bar and the inner bar, 
since the second mass inflow rate increases with the time intervals which are roughly the 
figure rotation period of the inner bar measured in the rotational frame of the outer bar. 

The high gas mass concentration cases in the small inner bars models are consistent 
with observations in our galaxy. In these models, a nuclear gas disk forms. Its size and its 
mass are R < 15 pc and ~ 1O 7 M , respectively. Interestingly, the size of the nuclear gas 
disk is very close to the location of the Arches cluster and the Quintuplet cluster. Moreover, 
the nuclear gas disk is massive enough to form these star clusters (we discuss this point in 
Sect J4.2p . Kinematics of gas induced by the inner bar is consistent with the molecular gas 
observations (we discuss this point in Sect l4.3p . On the other hand, in the small inner bar 
models with Qt < 0.05, the inner bar does not highly enhance mass inflow to the galactic 
center. Hence, the inner bar in our galaxy is Qt > 0.05, if mass supply to the Galactic 
center is due to the inner bar. 

There is difference betwe e n the size o f t he small inner bar models an d the inner bar 
reported by iNishiyama et al. tood . l2006h . INishiyama et al.l (120051 . 120061 1 trace ridge of 



distribution of red clump stars but do not show profile of gravitational potential of the inner 
bar. Our numerical results are consistent with their report, if non-axisymmetric component 
of gravitational potential of the inner bar is small beyond R = 200 pc. 



The large inner bars in our models are not consistent with observations in our galaxy, 
if mass supply to the galactic center is caused by the large inner bar. In the large inner 
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bar models with Qt > 0.11, high gas mass concentration occurs and the 200 pc gas ring is 
destroyed. This does not correspond to our galaxy. In the models with Qx < 0.11, the inner 
bar does not induce a large mass inflow to the galactic center. From these results, large inner 
bar is difficult to be the case in the Galaxy, if mass supply to the Galactic center is due to 
an inner bar. 

It is observed that velocity disper sion of gas clouds in the centra l region of the Galaxy is 



highe r than that in the Galactic disk (IRohlfs &: KreitschmanrJ 119871 ). lEnglmaier &: Gerhard 



19971 show that the gas flow can change drastically when the sound speed is changed, since 
existence and strengths of shocks depend on c s . In order to confirm the effect of the sound 
speed on the mass inflow, we try a test calculation in which the inner bar parameters are 
the same as the model S33 (f2iB = 300 km s _1 kpc" 1 ) and the artificial radial profile of c s , 
which rises from « lOkms" 1 at the inner edge of the 200 pc gas ring to 20kms~ 1 at the 
center, is assumed. The straight shocks in the inner bar become weaker and the mass inflow 
rate becomes smaller. We will study the effect of the sound speed on the gas flow further in 
a future work considering realistic cooling and heating process. 



4.2. Evolution of the nuclear gas disk 

We have shown that small massive gas disks form in the small inner bar models for > 
0.05 and their size are ~ 15 pc. It is interesting to study the self-gravitational instability of 
the nuclear gas disk. In an axisymmetric uniform thin gas disk, the dispersion relation of 
the small radial density perturbation in the axisymmetric mode is 

uj 2 = c 2 s k 2 - 2nGZ\k\ + k 2 , (3) 

where uj is the frequency of the perturbation, c s is the sound speed of gas, k is the wave 
number of the perturbation, G is the gravi tational constant, S is the surface density of the 



thin disk, and k is the epicyclic frequency ( iBinney &: Tremainelll987l ). From the dispersion 
relation, the density perturbation can grow if 

Q = -% < 1, (4) 
where Q is the Toomre Q- value. We define S crit as the surface density for Q — 1, 

c s (R)k(R) 

Scrit(-K) = -p, , (5) 

which may be the minimum surface density for the gravitational instability. Using S crit , we 
define M crit as 

r R 

M clit (R)= / 27rRE clit (R)dR' . (6) 
Jo 
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M cr it (R) may be a measure of gravitational instability of the disk. In the central several ten s 



parsecs of the galaxy, there is evidence for strong magnetic fields (e.g. IChuss et al.ll2003l ). 
Magnetic fields have an important role in the gravitational stability of the disk. To consider 
effect of the magnetic fields in the linear analysis, we assume simple configuration of the 
magnetic fields, since it is observationally unclear. We assume that the magnetic fields are 
parallel to the disk and homogeneous, B = Bpe ^, where Bp is a st rength of the magnetic 



fields and is the base vector of the azimuth. iFan fc Loul (119971 ) derived the dispersion 
relation 

u 2 = (c 2 s + c 2 A )k 2 - 2nGZ\k\ + k 2 (7) 
for this configuration, where c A is a Alfven velocity, 



CA 



B 2 
4iip 



We use this dispersion relation for our analysis. We assume that gas clumps are formed from 
perturbations with the largest growth rate. The wave length of the density perturbation 
with the largest growth rate is given by 

Amax - 2?F - 2C eff = 27TC ^(R) (g) 

where c e g = \fc^ + c 2 A . Gas clump mass is estimated as 

A m ax\ ^ 7I " 2c cfr(-^) 



M clump = vr(^) S crit = ^^. (10) 



Application this results to our numerical results shows that strong magnetic fields, which 
is comparable with the strongest magnetic fields observed in the Galactic central region, 
enable massive gas clumps to grow and these are comparable to the mass of the young 
massive star clusters in the Galactic center. Figure [10] show the result of the application for 
the model S33 (fijB = 300 km s" 1 kpc -1 ), which is one of the high gas mass concentration 
cases in the small inner bar models. In this figure, we assume that the gas in the nuclear 
gas disk sufficiently cools down to T = 100 K (c s ~ 1 km s _1 ). From this figure, the nuclear 
gas disk becomes gravitationally unstable after t = 300 Myr, if effect of the magnetic fields 
is very weak. The mass of the disk is 6.7 x 10 5 M Q at that time. The mass of the gas 
clumps is 100-300 M from equation (10). If Bq = 1 mG, the disk becomes unstable after 
t = 450 Myr. The mass of the disk is 2.9 x 10 6 M at that time. The mass of the gas clumps 
is 1.0-3.0 x 10 4 M . This mass is comparable to that of the young massive star clusters in 
the Galactic center. 
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Massive gas clumps can be formed even in the non-magnetic case. To investigate the 
non-linear evolution of the nuclear gas disk in the non-magnetic case, we perform very high 
resolution hydrodynamical simulations in paper II. In paper II, we show that many massive 
compact gas clumps are formed by gravitational instability of the cooling gas disk in the 
non-magnetic case. Typical mass and size of the clumps are several 10 3 M Q to 10 4 M and 
less than a few parsecs, respectively. The largest gas clumps have a mass of ~ 10 5 M . This 
is much larger than 100-300 M . This is because small gas clumps, which are formed rapidly 
from growth of density perturbation in the cooling disk, collide each other and merge into 
more massive clumps. The Arches and Quintuplet clusters have a mass of ~ 10 4 M and a 
size of < 1 pc. If we assume a star formation efficiency of ~ 0.1, these clusters can be formed 
from the gas clump of mass ~ 10 5 M , which is comparable to the largest gas clumps in our 
numerical results of paper II. 



4.3. Longitude-velocity diagrams of gas flow in the nested bars 

We make longitude- velocity (l-v) diagrams from our numerical results with the following 
two aims. One is to compare our numerical results with observations in our galaxy. Another 
is to show that characteristic features of gas motion induced by the inner bar can be evidence 
of inner bars in external galaxies. 

Figure [TT1 shows the l-v diagrams of the model S33 (ilm = 300 km s" 1 kpc" 1 ), which is 
one of the high gas mass concentration cases in the small inner bar models, for Ad' = 0° 
and 90°, where A9' is the angle between the direction of the inner bar and the Sun-Galactic 
center line. In the diagram, we assume that the outer bar is inclined at an angle of 20° with 
respect to the Sun-Galactic center line, that the distance of the Sun from the Galactic center 
is 8kpc, and th at the circular velocity of the Sun is 220 km s -1 . These assumption is based 



on the results of iBissantz et al.l (120031 ). In this figure, we classify the gas components of the 
results into 7 groups by colors (the detail of the classification is described in the caption of 
Fig. [TT]) according to property of gas motion. The nuclear gas disk component is shown by 
the red points in Fig. [TTJ The nuclear gas disk component in Fig. [11] is weakly dependent on 
AO', since circular motions dominate in the disk. The straight shocks component is shown 
by the green points in Fig. [TTJ The feature of the straight shocks depends on AO'. When 
the inner bar is perpendicular to the outer bar, the straight shocks component is clearly 
distinguishable from the nuclear gas disk component and the 200 pc gas ring component. 
The elliptical gas ring component is shown by the purple points in Fig. [TTJ The feature of 
this ring strongly depends on A8', since it is elongated along the inner bar. 



There are many observational studies on gas distribution and kinematics in the central 
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region of our galaxy. IStark et all (120041 ) give the l-v diagram of highly excited rotational 
emission lines of CO (J=4-3 and J=7-6) in the central region of our galaxy observed by 
AST/RO. They cover a range of —1.2° < I < 2°. Their l -v dia gram traces high density 
components of molecular gas. iRodriguez-Fernandez et all (120061 ) show the l-v diagram of 
CO (J=2-l ) using the publish ed data. Their l-v diagram covers the same region as the l-v 
diagram of IStark et all (120041 ). but traces diffuse molecular gas. It is known that there are 
two compact GMCs, the 20 km s" 1 cloud and the 50 km s~ x cloud in the Galactic center 
region. The 20kms _1 cloud is located at R < 10 pc from the center in the projection. It 
has a total mass of ~ 3 x 10 5 M and its radi al velocities is in a range of ~ 5-25 km s" 1 . 
The 50 km s _1 cloud have a mass of ~ 10 5 (IMezger et al.lll996h. The position s of these 
GMCs in the l-v diagram are shown in iNagayama et all J2007h ~ lOka et al.l (120071 ) give the 
l-v diagram of a highly excited rotational emission line of CO (J=3-2) with high resolution 
from / = —0.2° to 0.1°. They show that there is a pair of high velocity emission (they are 
named CND + and CND~ in their paper) within 0.05° ~ 6.5 pc from Sgr A*. The line-of-sight 
velocity of CND + and CND~ is 50-100 km s -1 and —50 — 120 km s" 1 , respectively. 

Our numerical results is consistent with the observations in the central region of our 
galaxy. The nuclear gas disk component in Fig. [TT] for A6' = 90° is in the longitude range 
of —0.2° < / < 0.2° and in the velocity rang e of —100 km s -1 < v < 100 km s -1 . This is 
same range of the most inner orbit shown in I St ark et all (120041 ) . The velocity range of the 
nuclear gas disk agrees with that of the CND. The similar agree ment between the nuclear gas 



(2006 


) and 


Oka et al. 


(2007) 



and 50 km s 1 clo uds lie in the same region as 



the nuclear gas disk in the l-v diagram (INagayama et all 120071 1 . Thus, the nuclear gas disk 
component well corresponds to the observations. There are not clear high velocity compo- 
ne nts corresponding to th e elliptical gas ring component fo r Ad' = 0° in the l-v diagrams 
of IStark et al.l (120041 ) and iRodriguez- Fernandez et al.l (120061 ). The other gas components in 
our l-v diagrams occupy the same region in their l-v diagrams. Thus, our l-v diagrams for 
AO' = 90° well corresponds to our Galaxy. 

We compare our numerical results with molecular gas observation in Maffei 2. Our l-v 
diagram at A9' = 90° w ell corresponds to t he CO position-ve locity (p-v) diagrams of the 
nuclear region of Maffei 2 (IMeier et al.ll2008l ). iMeier et all (120081 ) performed an observation of 
the nuclear region of Maffei 2 with high spacial resolution with the OVRO and BIMA arrays 
and found a parallelogram feature and two intense features at both side of the parallelogram 
feature in their p-v diagrams. The parallelogram feature extends over —5" < p < 15" and 
—125 km s" 1 < v < 125 km s" 1 in their diagrams. The two intense features are located 
at (p,v) ~ (—15", 50 km s -1 ) and (20", —50 km s" 1 ) in their diagram . They explain these 
features by simple linear orbits in their nuclear bar model. The nuclear gas disk component 
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in our l-v diagram corresponds to the parallelogram feature. The straight shocks component 
and the elliptical gas component in our l-v diagram well correspond to the two intense 
features. Thus, our results strongly support their interpretation that Maffei 2 likely has an 
nuclear bar. We propose that the nuclear gas disk component, the straight shocks component, 
and the elliptical gas component are indirect evidence for an inner bar. Observation of 
molecular gas in the nuclear region of external barred galaxies with high spacial resolution, 
e.g. ALMA, can give evidence of inner bars, even if they are hidden by a large amount of 
gas and dust. 



4.4. Important role of central mass concentration 



We discuss the difference between our numerical results and numerical results of lMaciejewski et al 



(120021 ). In our simulations, the massive nuclear gas disks are formed in the galactic cen- 
ter. Formation of the nuclear gas disks is due to the straight shocks inside the inner bars. 
On the other hand, bot h such nuclear gas disks and straight shocks are not formed in 
Maciejewski et al.l ( 120021 ). although they also simulated gas flow in a nested barred model. 



We consider the central mass concentration as the main reason for the difference, since 
major difference between our models and their model is central mass concentration. We 
assumed the high cen tral mass concen t ration that is modeled on the basis of the nuclear 
bulge pr ofile given by lLaunhardt et al.l (120021 ). while th e central mass concentration in the 
model of lMaciejewski et al.l (120021 ) is low (see Fig. 3 in iMaciejewski fc Sparkdl2000l ). It is 
shown that a high central mass concentration in a ba r red potential strongly affects orbital 
structure of stars and gas (IFukuda et al.l Il998l . l2000l ; lAnn &: Thakurl 120051 ) . The central 
mass concentration tends to change the shape of the orbits of stars into rounder shapes at 
the nearer central region of the galaxy. When the galaxy has an inner bar, the shapes of 
the orbits are elongated at the radii which are comparable to the semi-major axis of the 
inner bar. In smaller radii, the shape of the orbits changes into the circular orbits in the 
inner bar potential, if the central mass concentration is sufficiently high. Straight shocks 
may form, if the shape of the orbits rapidly vary as the radii becomes small, since gas 
collides each other at the region where the orbits are overcrowded and dissipates. Hence, a 
high central mass concentration is important for formation of straight shocks and therefore 
formation of nuclear gas disks. We conclude that the difference between our numerical results 
and numerical results of lMaciejewski et al.l (120021 ) is mainly due to the difference in a central 
mass concentration. It is important to study a self-consistent model of nested barred galaxies 
with high central mass concentrations and their stability. 
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5. SUMMARY 

We summarize our study as follows: 

. We have performed two dimensional hydrodynamical simulations to investigate mass 
supply process by nested bars. We have a ssumed the gravitatio nal potential model of 



our galaxy, base d on the Galaxy models of iBissantz et all (120031 ) and the nuclear bulge 



profile given by lLaunhardt et al.l (120021 ) adding an inner bar. We have assumed two 



cases of the size of the inner bar models, oib = 200 pc and 600 pc. 

2. In the small inner bar models, a large amount of gas concentrates into the galactic 
center for 1) 0.05 < Q T < 0.12 and tt m ~ 225 km s^ 1 kpc" 1 and 2) Q T > 0.12 and, 
JIib ~ (f2 — ft/2) max or Qib ~ 225 km s _1 kpc -1 . The straight shocks are formed within 
the inner bar. This is partly due to that Qt in these models is high and partly due 
to that the central mass concentration in our models is high. The straight shocks 
sweep gas in the inner bar region. The gas trapped by the straight shocks falls into 
the galactic center and then the nuclear gas disk is formed at the center. The size and 
mass of the nuclear gas disk are < 15 pc and ~ 10 7 M , respectively. 

3. In the large inner bar models, a large amount of gas concentrates into the galactic 
center for Qt > 0.11. In the course of the gas concentration, the inner bar destroys 
the 200 pc gas ring. The destruction of the 200 pc gas ring is not consistent with the 
CMZ. We conclude that the inner bar of our galaxy is not both large and strong, if 
recent mass supply to the galactic center is due to the inner bar of our galaxy. 

4. The high gas mass concentration cases in the small inner bar models well agree with 
the observed feature as follows. Extent and kinematics of the nuclear gas disk in our 
results are consistent with the observations of the molecular gas in the central region of 
our galaxy. The size of the nuclear gas disk is very close to the location of the Arches 
cluster and the Quintuplet cluster, and its mass is enough to form these star clusters. 

5. We have discussed the self-gravitational instability of the nuclear gas disk formed in 
our simulations. Assuming magnetic fields as strong as observed one in the central 
tens parsecs of our galaxy, the most rapid growing unstable mode corresponds to gas 
clumps which have comparable mass to the Arches and Quintuplet cluster. In next 
paper, we will study non-linear evolution of massive nuclear gas disks. 

6. We have shown the characteristic features in the l-v diagram induced by the small 
inner bar. These features can be clues about existence of inner bars in extra galaxies. 
They will be useful for future observation of central regions of galaxies, e.g. ALMA. 
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Fig. 1. — Rotation curves in the entire (left panel) and the central region (right panel) of 
the model S33. These lines show rotation curves by total mass (solid line), disk and bulge 
(dashed line), inner bar (dotted line), and dark halo (dotted- dashed line), respectively. 

The high resolutional version of the figure is available from 
http : / /astro3 . sci . hokudai . ac . jp/~name/. 



This preprint was prepared with the AAS IATgX macros v5.0. 
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Fig. 2. — Angular frequency curve in the entire (left) and the central region (right). These 
lines show angular frequency curves of Q — nj2(solid line), fl(dashed line), and Q + hi/2(dotted 
line) and the pattern speed of the outer bar QoB(dash- dotted line). The positions of the 
inner Lindblad resonance (-Rilr), the corotation resonance (Rcr), and the outer Lindblad 
resonance (-Rolr) of the outer bar are 750 pc, 3750 pc and 6750 pc, respectively. 

The high resolutional version of the figure is available from 
http : / /astro3 . sci . hokudai . ac . jp/~name/. 



Table 1. Axial ratio and mass of the small inner bar models. 
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(0.218) 
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S35 
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(0.274) 
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S24 
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(0.036) 


(0.048) 


(0.071) 


(0.115) 


4/3 






S13 


S14 


S15 








(0.012) 


(0.019) 


(0.031) 



a The mass of the inner bar model. 

b Thc axial ratio of the inner bar model, aie is a semi-major axis of the 
inner bar model and 6ib is a semi-minor axis of the inner bar model. 

c The values in parentheses show Qt- 
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Fig. 3. — Surface density of gas in the model N at t — 500 Myr. The gray-scale bar show the 
surface density of the range of 1 M Q pc _2 -10 5 M Q pc -2 in the logarithmic scale. The outer 
bar lies in the horizontal direction in all the panels. 

The high resolutional version of the figure is available from 
http : / /astro3 . sci . hokudai . ac . jp/~name/. 



Table 2. Axial ratio and mass of the large inner bar models. 









M IB [M Q ] a 






aiB/6lB b 


1.0 x 10 8 


2.5 x 10 8 
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1.0 x 10 9 
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(0.276) 


2 
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(0.036) 


L23 
(0.071) 


L24 
(0.105) 


L25 
(0.138) 


4/3 






L13 
(0.023) 


L14 
(0.036) 


L15 
(0.046) 



a The mass of the inner bar model. 

b The axial ratio of the inner bar model, aie is a semi-major axis of the 
inner bar model and 6jb is a semi-minor axis of the inner bar model. 

c The values in parentheses show Qt- 
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Fig. 4. — Time variation of the surface density of the central kpc of the model S33 
(300 km s" 1 kpc -1 ). The gray-scale range is the same as Fig. [3j We show the calcula- 
tion time and the angle between the outer bar and the inner bar, A6, in the upper right 
corner in the each panel. 



The high resolutional version of the figure is available from 
http : / /astro3 . sci . hokudai . ac . jp/~name/. 
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Fig. 5. — Velocity fields of the model S33 (300 km s 1 kpc 1 ). The gray-scale are the same 
as FigfJ Upper panels : A0 = 0° (t = 490 Myr). Lower panels : A9 = 90° (t = 497 Myr). 
Note that the arrows show velocities in these region in relative scale. 



The high resolutional version of the figure is available from 
http : / /astro3 . sci . hokudai . ac . jp/~name/. 
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Fig. 6. — Surface density distribution of the model S21. Upper row : the model S21 {Vim = 
300 km s -1 kpc" 1 ). Lower row : the model S21 (Qm = 200 km s _1 kpc" 1 ). First and second 
columns : AO = 0°. Third and fourth columns : A6 = 90°. 

The high resolutional version of the figure is available from 
http : / /astro3 . sci . hokudai . ac . jp/~name/. 
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Fig. 7.— M 20 (t) of different three models, S33 (n 1B = 300 km s" 1 kpc" 1 ), S21 (Q m = 
300 km s -1 kpc" 1 ), and N (upper left panel). M 2 o(t) of the model S21 (upper right panel) 
and the model S33 (lower panel) for various pattern speeds of the inner bar. The horizontal 
and vertical axis are calculation time and logarithmic value of M2o(t)/M & , respectively. 



The high resolutional version of the figure is available from 
http : / /astro3 . sci . hokudai . ac . jp/~name/. 
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Fig. 8. — Time variation of the surface density of gas in the model L42. The gray-scale bar 
shows the surface density of the range 1 M pc _2 -10 6 M pc -2 in the logarithmic scale. The 
values of upper right corner in each panel are the same as Fig]H 

The high resolutional version of the figure is available from 
http : / /astro3 . sci . hokudai . ac . jp/~name/. 
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Fig. 9. — Gas mass concentration in our numerical results of the small inner bar models. The 
horizontal and vertical axis are Qt and Qib, respectively. The high gas mass concentration 
models are shown by the filled circles, the low gas mass concentration models are shown by 
the crosses. 



The high resolutional version of the figure is available from 
http : / /astro3 . sci . hokudai . ac . jp/~name/. 



Table 3. Pattern speeds of the inner bar models. 
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a The open circles show the model we simulate. 
b The letter 'L' means the large inner bar models. 



Table 4. High gas mass concentration cases in the small inner bar models. 



Pattern speed of the inner bar Htb [km s 1 kpc "*"] 
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l"The models in which the second mass inflow starts just before the end of the simulation. 
*Thc models in which the second mass flow rate changes periodically. 



Table 5. High gas mass concentration cases in the large inner bar models. 



Model Name 
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Fig. 10.— Left: Time evolution of M(< R) of the model S33 (300 km s" 1 kpc" 1 ) and the 
critical disk mass M cr i t (R). M(< R) at each time are shown by dashed lines, M c[Tt (R) are 
shown by two solid lines for the case of B = mG and B = 1 mG. Right: The mass of the 
gas clump M c i ump given by eq.(10) at each radius. 



The high resolutional version of the figure is available from 
http : / /astro3 . sci . hokudai . ac . jp/~name/. 



- 36 - 



250 
200 
150 
100 
50 
u 
-50 
-100 
-150 
-200 
-250 

















\ 


■ 








"' ■■■ 








\ 






i , . 




2 1 D -J 


2 




I - ll v. ' I: 1 [ii'il'l-k !-!*_=.- 1 -->_ 1 





250 
200 
150 
100 

-~ 50 
g 

-± 

I "5" 
-ICO 
-L50 
-200 
-250 



^-;. : ' : :?L^3;!fa»^ ... _ 








m 
















































• 













l n -i 



Fig. 11. — Longitude-velocity diagram of the model S33 (£lm = 300 km s -1 kpc" 1 ). The 
horizontal and vertical axis represent Galactic longitude [°] and line-of-sight velocity [km s -1 ], 
respectively. Left panel shows the l-v diagram for A6' = 0° (t = 489 Myr). Right panel shows 
the l-v diagram for A6' = 90° (t = 495 Myr). Red points: the nuclear gas disk component; 
green points: the straight shocks component; blue points: the inner void component, which 
corresponds to low density regions within the elliptical gas ring; purple points: the elliptical 
gas ring component; aqua points: the outer void component, which corresponds to low 
density regions between the elliptical ring and the 200 pc gas ring; yellow points: the 200 pc 
gas ring component; black points: the region of R > 300 pc. 



The high resolutional version of the figure is available from 
http : / /astro3 . sci . hokudai . ac . jp/~name/. 



